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Prior  Merging  Techniques 


Numerous  Previous  Merge  Methods: 
o  2:1  -  Cannot  Conserve  Energy 

(Lapenta  &  Brackbill,  JCP  1994) 

o  Complex  Macro-particles  with 
Internal  Energy 

(Hewett,  JCP  2003) 
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OD-Merge  Examples 
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Increased  Complexity  First  Attempt... 

o  Annular  Potential: 

$  oc  log(r / r-m)  / log(rm/ rout) 
o  Non-Trivial  Cartesian  Derivative  of 
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Annular  Test  Case 


Increased  Complexity  First  Attempt... 


9  Annular  Potential: 

$  oc  log(r/rm)/log(riB/r0U t) 

9  Non-Trivial  Cartesian  Derivative  of  <f> 

9  Stable  Spiral  Required  C.N.  Push 

9  Different  Periods  per  Particle  - 
Non-Periodic  unlike  Parabolic  Well 

o  Part./Cell  Drops  &  #-Cells  Active  Grows 
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Annular  Test  Case 


Increased  Complexity  First  Attempt... 


9  Annular  Potential: 

$  oc  log(r / r-m) / log(rm/ rout) 

9  Non-Trivial  Cartesian  Derivative  of  <f> 

9  Stable  Spiral  Required  C.N.  Push 

9  Different  Periods  per  Particle  - 
Non-Periodic  unlike  Parabolic  Well 

o  Part./Cell  Drops  &  #-Cells  Active  Grows 

Not  a  Good  Repeating  Merge/Split  Test 
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Spherical  Test  Case 


Increased  Complexity  Second  Attempt... 


o  Spherical  Potential:  &  ocl/r 
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Spherical  Test  Case 


Increased  Complexity  Second  Attempt... 

o  Spherical  Potential:  oc  1  jr 
©  Slight  AX  — »  Orbits  Unsynchronized 
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Without  Merge  &  Split 
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Spherical  Test  Case 
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Spherical  Test  Case 


Increased  Complexity  Second  Attempt... 

o  Spherical  Potential:  oc  1  jr 
o  Slight  AX  Orbits  Unsynchronized 
©  Periods  Synched  via  Orbital  Mech.  Eqns. 
©  Periodic  with  Full  Nonlinear-C.N.  Push 
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Without  Merge  &  Split 


■ 


Spherical  Test  Case 


Increased  Complexity  Second  Attempt... 

o  Spherical  Potential:  <1>  oc  1/r 
o  Slight  AX  — )►  Orbits  Unsynchronized 
©  Periods  Synched  via  Orbital  Mech.  Eqns. 
©  Periodic  with  Full  Nonlinear-C.N.  Push 

©  Adding  Merge  &  Split: 

Energy  Conservation  OK 
Accuracy:  Scatter  +  Dispersion? 


Including  Merge  &  Split 
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Spherical  Test  Case 


Increased  Complexity  Second  Attempt... 

o  Spherical  Potential:  <f>  oc  1/r 
o  Slight  AX  — »  Orbits  Unsynchronized 
©  Periods  Synched  via  Orbital  Mech.  Eqns. 
©  Periodic  with  Full  Nonlinear-C.N.  Push 

©  Adding  Merge  &  Split: 

Energy  Conservation  OK 
Accuracy:  Scatter  +  Dispersion? 

Source  of  Dispersion? 
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Dispersion  from  Incomplete  2nd  Moment?  i 


0(Ax^ell)  Potential  Conservation: 

o  Remap  must  Preserve  Identical 
Contraction,  8$$  :  AxtAxj 


'®xx 

®Xy 

3>xz~ 

~AxAx 

AxAy 

AxAz 

&yx 

<f>yy 

3>yz 

AyAx 

AyAy 

AyAz 

*zy 

®zz\ 

_AzAx 

AzAy 

AzAz_ 

d 


□  ►  <  g  >  i  =  ►  <  ^  >  jf  c\  o 


R.S.  Martin  (ERC  Inc.) 


Distribution  A:  Approved  for  public  release;  unlimited  distribution;  PA  #TBD  14  / 19 


Dispersion  from 


0(Ax^ell)  Potential  Conservation: 

o  Remap  must  Preserve  Identical 
Contraction,  8$$  :  AxtA. xj 

©  2-Particle  Only  Reproduces  Diagonal 
with  x\a/b^  =x~i±ei  if  a  =  Axf/2 
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Dispersion  from 


0(Ax^ell)  Potential  Conservation: 

o  Remap  must  Preserve  Identical 
Contraction,  8$$  :  AxtA. xj 

©  2-Particle  Only  Reproduces  Diagonal 
with  x\a/b^  =Xi±€i  if  a  =  Axf/2 

©  Non-Degenerate  aXiaXj  -A  3 -Eigenvalues 
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2nd  Moment?  ijfe 


'Ax  Ax  Ax  Ay 
AjAx  Ay  Ay 
_AzAx  AzAy 


AxA  z 
AyAz 
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Dispersion  from  Incomplete  2nd  Moment?  o 


0(Ax^ell)  Potential  Conservation: 

o  Remap  must  Preserve  Identical 
Contraction,  8$$  :  AxtA. xj 

©  2-Particle  Only  Reproduces  Diagonal 
with  x\a/b^  =Xi±€i  if  a  =  Axf/2 

©  Non-Degenerate  aXiaXj  -A  3 -Eigenvalues 
©  Non-Deg.  — >•  4+  Non-Coplanar  Particles 

(3-Particles  =  Plane  &  A3  =  0;  2-Particles  =  Line  &  A2 ,  A3  =  0) 


'Ax  Ax  Ax  Ay  AxAz' 
AyAx  Ay  Ay  AyAz 
_AzAx  AzAy  AzAz_ 
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Dispersion  from  Incomplete  2nd  Moment?  o 


0(Ax^ell)  Potential  Conservation: 

9  Remap  must  Preserve  Identical 
Contraction,  8$$  :  AxtA. xj 

9  2-Particle  Only  Reproduces  Diagonal 
with  x\a/b^  =Xi±€i  if  a  =  Axf/2 

9  Non-Degenerate  aXiaXj  -A  3 -Eigenvalues 
9  Non-Deg.  — >•  4+  Non-Coplanar  Particles 

(3-Particles  =  Plane  &  A3  =  0;  2-Particles  =  Line  &  A2 ,  A3  =  0) 

9  Same  — >«  4  Merge  for  2nd -Vq\.  Moment 
9  “Principal  Component  Analysis” 


'Ax  Ax  Ax  Ay  AxAz' 
AyAx  Ay  Ay  AyAz 
_AzAx  AzAy  AzAz_ 
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Dispersion  from  Incomplete  2nd  Moment?  o 


0(Ax^ell)  Potential  Conservation: 

o  Remap  must  Preserve  Identical 
Contraction,  :  AxtA xj 

9  2-Particle  Only  Reproduces  Diagonal 
with  x\a/b^  =Xi±ei  if  €i  = 

9  Non-Degenerate  aXiaXj  -A  3 -Eigenvalues 
9  Non-Deg.  — >■  4+  Non-Coplanar  Particles 

(3-Particles  =  Plane  &  A3  =  0;  2-Particles  =  Line  &  A2 ,  A3  =  0) 

9  Same  4  Merge  for  2nd -Vq\.  Moment 

9  “Principal  Component  Analysis” 

9  Improved  but  Non-Negligible  Scatter  and 
Dispersion 
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Dispersion  from  XV  Scatter? 
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Kinetic  Potential  Energy  Scatter: 
o  A  Closer  look  at  Merge  — >•  2 
o  «  Match  sgn(efv))  to  sgn(AviAvj)  Moments 
o  Select  sgn^x)^  to  Match  s^a^Av,) 
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Dispersion  from  XV  Scatter? 


Kinetic  Potential  Energy  Scatter: 


o  A  Closer  look  at  Merge  — >•  2 
o  «  Match  sgn(ejv))  to  sgn(AViAvj)  Moments 
O  Select  Sgn(ef})  to  Match  Sgn(AxjAvj) 

©  Some  Improvement  Scatter  &  Dispersion 
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Dispersion  from  XV  Scatter? 


0 


Kinetic  Potential  Energy  Scatter: 


9  A  Closer  look  at  Merge  — >•  2 
9  w  Match  sgn(e^v))  to  sgn(AViAvj)  Moments 
9  Select  sgnfaW)  to  Match  sgn(AXiAvi) 

9  Some  Improvement  Scatter  &  Dispersion 


9  Why  Stop  at  just  XV  Sign? 

=  (A^)  /eW 


4 
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Kinetic  Potential  Energy  Scatter: 

©  A  Closer  look  at  Merge  — >•  2 
©  «  Match  sgn(efv))  to  sgn(AviAvj)  Moments 
©  Select  Sgn(e^  to  Match  Sgn(AXiAvi) 

©  Some  Improvement  Scatter  &  Dispersion 
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©  Why  Stop  at  just  XV  Sign? 


=  (A^)  /Ci(v) 


©  Much  Improved  Scatter  &  Dispersion 
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Dispersion  from  XV  Scatter? 


Kinetic  Potential  Energy  Scatter: 

Q  A  Closer  look  at  Merge  — >•  2 
o  «  Match  s8n^v)^  to  sgn(AviAvj)  Moments 
o  Select  sgn^x)^  to  Match  ^(Ax^AvO 
o  Some  Improvement  Scatter  &  Dispersion 
o  Why  Stop  at  just  XV  Sign? 

6«=(S^)/CM 

©  Much  Improved  Scatter  &  Dispersion 
©  Sphere  Sensitive:  Kinetic  <->  Potential 


Dispersion  from  XV  Scatter? 


Kinetic  Potential  Energy  Scatter: 


o  A  Closer  look  at  Merge  — >•  2 
o  w  Match  sgn(e^v))  to  sgn(AViAvj)  Moments 
O  Select  Sgn(efx))  to  Match  Sgn(AXiAVi) 

©  Some  Improvement  Scatter  &  Dispersion 

©  Why  Stop  at  just  XV  Sign? 

=  (A^)  /Ci(v) 

©  Much  Improved  Scatter  &  Dispersion 

©  Sphere  Sensitive:  Kinetic  <->  Potential 

©  K.E.  &  P.E.  Conserved  Independently... 
XV-Moment  Better  Combination 
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Review  Sphere  Test  Results 
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Merge  Results  after  1 -Cycle 

•  ••  fj] '  v  '  t 

No  Merge  Random  Sign  PC  A  XV-Sign  XV-Full 


Note:  Results  Due  to  Extreme  Case  Sensitivity 
C.N.  Push  Required  for  Stability 
Period  Sync  Required  for  Initial  Conditions 


Review  Sphere  Test  Results 
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Merge  Results  after  1 -Cycle 


Q>  'y  (J>  «»* 

Not  Synched  Random  Sign  PCA  XV-Sign  XV-Full 


Note:  Results  Due  to  Extreme  Case  Sensitivity 
C.N.  Push  Required  for  Stability 
Period  Sync  Required  for  Initial  Conditions 


Review  Sphere  Test  Results 


Merge  Results  after  1 -Cycle 

Q  j  ‘F  <J>  '»* 

Not  Synched  Random  Sign  PCA  XV-Sign  XV-Full 

Note:  Results  Due  to  Extreme  Case  Sensitivity 
C.N.  Push  Required  for  Stability 
Period  Sync  Required  for  Initial  Conditions 
Result  of  Minor  Velocity  Difference  Across  IC  Box 


Review  Sphere  Test  Results 


0 


Merge  Results  after  1 -Cycle 

Q  j  ‘F  <J>  '»* 

Not  Synched  Random  Sign  PCA  XV-Sign  XV-Full 

Note:  Results  Due  to  Extreme  Case  Sensitivity 
C.N.  Push  Required  for  Stability 
Period  Sync  Required  for  Initial  Conditions 
Result  of  Minor  Velocity  Difference  Across  IC  Box 
Merge  Scatter  in  Cells  <  Synch  A  V 
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Future  Work:  Combined  PCA /XV  Merge 
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Future  Work:  Combined  PCA /XV  Merge 


PCA  +  XV-Merge? 

o  Both  Partially  Inhibited  Dispersion 
®  Extra  Degrees  of  Freedom  in  PCA 
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Future  Work:  Combd* 


PCA  +  XV-Merge? 

o  Both  Partially  Inhibited  Dispersion 
©  Extra  Degrees  of  Freedom  in  PCA 
«  Use  Extra  DOF  for  XV-Match? 
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Future  Work:  Combined  PCA /XV  Merge 


PCA  +  XV-Merge? 

o  Both  Partially  Inhibited  Dispersion 

o  Extra  Degrees  of  Freedom  in  PCA 

«  Use  Extra  DOF  for  XV-Match? 

o  XV-Dispersion  Source  of  Error  in 
Lapenta’s  MHD  Slow  Shock? 


Lapenta  -  JCP  181,  317-337  (2002) 
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AFRL/RQRS  M&S  Future  Work 


0 


Integrate  R&D  w/  Production  TODO: 
o  Combine  PCA  and  XV  Merges 

(Martin/Lederman) 

o  Remap  <->  Vlasov  <->  Fluid 

(Martin/Bilyeu/TBD) 

o  Implicit/Multiscale  GPU-Accel.  PIC 

(Lederman/Gimelsheins/Martin/TBD) 

o  GPU  Accelerated  Chemical  Kinetics  / 
CR  Ar-Ne-Xe-Molecular  Models 

(Le/Cole*/Kapper*  PhD  Research) 


*Note:  Former  Co-op  Student  Work  to  be  Integrated  into  Framework 
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